Phase detection at the quantum limit with multi-photon Mach-Zehnder interferometry 
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We study a Mach-Zehnder interferometer fed by a coherent state in one input port and vacuum 
in the other. We explore a Bayesian phase estimation strategy to demonstrate that it is possible 
to achieve the standard quantum limit independently from the true value of the phase shift and 
specific assumptions on the noise of the interferometer. We have been able to implement the protocol 
using parallel operation of two photon-number-resolving detectors and multiphoton coincidence logic 
electronics at the output ports of a weakly- illuminated Mach-Zehnder interferometer. This protocol 
is unbiased and saturates the Cramer-Rao phase uncertainty bound and, therefore, is an optimal 
phase estimation strategy. 

PACS numbers: 42.50 St, 42.50 -p 



The Mach-Zehnder (MZ) interferometer [H,[l] is a truly 
ubiquitous device that has been implemented using pho- 
tons, electrons [3j, and atoms [4|, [5[. Its applications 
range from micro- to macro-scales, including models of 
aerodynamics structures, near-field scanning microscopy 
and the measurement of gravity accelerations Q • The 
central goal of interferometry is to estimate phases with 
the highest possible confidence [1, H, Q3] while taking into 
account sources of noise. Recent technological advances 
make it possible to reduce or compensate the classical 
noise to the level where a different and irreducible source 
of uncertainty becomes dominant: the quantum noise. 
Given a finite energy resource, quantum uncertainty prin- 
ciples and back reactions limit the ultimate precision of 
a phase measurement. In the standard configuration of 
the MZ interferometer, a coherent optical state with an 
average number of photons n = \a\ 2 enters input port 
a and the vacuum enters input port 6, as illustrated in 
Figure (pQ) . The goal is to estimate the value of the phase 
shift after measuring a certain number of photons N c 
and Nd at output ports c and <i, which, in the experi- 
ment discussed in this Letter, is made possible by two 
number-resolving photodetectors. 

The conventional phase inference pro tocol estimates 
the true value of the phase shift as [ll|, [TH, [l3[ : 
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difference detected at the output ports, averaged over p 
independent measurements. The phase uncertainty of 
estimator (pQ) is 
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which follows from a linear error propagation theory. 
Eq. ([2]) predicts an optimal working point at a phase shift 
6 = 7r/2, where the average photon number difference 
varies most quickly with phase. As 6 approaches or 7r, 
the confidence of the measurement becomes very low and 
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FIG. 1: (a) Schematic of a Mach-Zehnder interferometer. A 
phase sensitive measurement is provided by the detection of 
the number of particles N c and Nd at the two output ports, 
(b) Pulse height distribution for a visible light photon counter 
(VLPC) used in the experiment. The power incident on the 
detector is 144 fW at a wavelength of 780 nm. The vertical 
lines show the decision thresholds 231. 



eventually vanishes. As a consequence, this interferomet- 
ric protocol does not allow the measurement of arbitrary 
phase shifts. This can be a serious drawback for applica- 
tions like laser gyroscopes, the synchronization of clocks, 
or the alignment of reference frames. Furthermore, to es- 
timate small phase shifts with the highest resolution, the 
interferometer has to be actively stabilized around tt/2. 
This generally requires the addition of a feedback loop, 
which can be quite costly in terms of time and energy 
resources. 

It was first noticed by Yurke, McCall and Klauder 
(YMK) in [14[ that the estimator (pQ) does not take into 
account all the available information, and, in particular, 
the fluctuations in the total number of photons at the 
output ports. The possibility to improve Eq.([2|) is con- 
firmed by the analysis of the Cramer-Rao lower bound 
(CRLB) 0, [l6|], which provides, given an input state 
and choice of observables, the lowest uncertainty allowed 
by Quantum Mechanics. For a generic, unbiased, estima- 
tor, AQcrlb = l/y/pF(0), where F(0) is the Fisher in- 
formation [17j, which, in general, can depend on the true 
value of the phase shift 6 and the number of independent 
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FIG. 2: (color online). Phase distribution P((j)\N c , N d ), for: 
A) iV c = 0, N d = 0; B) N c = 0, N d = 1; C) iV c = 1, JV d = 1; 
D) iV c = 0, N d — 2. The circles are the experimental data 
collected in the calibration part of the experiment, the dashed 
line is the ideal phase distribution Eq.Q. The solid line is a 
fit of the data according to Eq.([5]). 



FIG. 3: (color online). Phase probability distribution Eq.©, 
obtained after p independent experimental measurements 
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1, B) p = 10, C) p = 100, D) 



p — 1000. The true value of the phase shift is 0/tt 
shown by the vertical dashed line. 
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measurements p. A direct calculation of the Fisher in- 
formation for the coherent (g) vacuum input state, gives 
F(6) = n. Therefore, the Cramer-Rao lower bound is 
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which, in contrast with the result of Eq.©, is indepen- 
dent of the true value of the phase shift. The only as- 
sumption here is that the observable measured at the out- 
put ports is the number of particles. It is well known (see 
for instance [18] ) that the Maximum Likelihood (ML) 
estimator, defined as the maximum, Oml, of the Like- 
lihood function P(iV c , Nd\4>) (see below), saturates the 
CRLB, but only asymptotically in the number of mea- 
surements p. In the current literature there have been 
alternative suggestions to obtain an unbiased estima- 
tor and a phase independent sensitivity with a Mach- 
Zehnder interferometer 0, Q H, H El • They will be 
specifically addressed at the end of this Letter. 

Here we develop a protocol based on a Bayesian 
analysis of the measurement results [1, [H, 
EH. The goal is to determine P((j)\ 7V C , N d ), 
the probability that the phase equals <p given 
the measured N c and N d . Bayes' theorem pro- 
videsthis: P(</>\N c ,N d ) = P(N c ,N d \</>) P(</>)/P(N c ,N d ), 
where P(N C: N d \<p) is the probability to detect N c and 
N d when the phase is (j) [22], P{4>) quantifies our prior 
knowledge about the true value of the phase shift, and 
P(N Cl N d ) is fixed by normalization. Assuming no prior 
knowledge of the phase shift, P(4>) = 1/tt. In the ideal 
case, the Bayesian phase probability distribution can be 
calculated analytically for any value of N c and N d , 



P(^|iV C) iV d )=c(cos^) 2JVc (sin^) 
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where C 



r " 1 (i + C Nc +n 2 (T) Nd " ^ s a normalization con- 
stant. In practice, one must measure P(N c ,N d \<j)) and, 
from this, determine P{(j)\N c , N d ). This distribution pro- 



vides both an estimate on the phase and the uncertainty 
in this estimate. 

There are several advantages to using a Bayesian pro- 
tocol. Notably, it can be applied to any number p of 
independent measurements, it does not require statisti- 
cal convergence or averaging, and it provides uncertainty 
estimates tailored to the specific measurement results. 
For instance, with a single measurement, p — 1, it pre- 



dicts an uncertainty that scales as A6 
Eq.(j4j) does not depend on n, 



l 



Since 



VN C +N d ' 

the estimation is insensitive 
to fluctuations of the input laser intensity. Most impor- 
tantly, its uncertainty, in the limit p >> 1, is AO = —^== 
which coincides with the CRLB. 

To implement the proposed protocol we have real- 
ized a polarization Mach-Zehnder interferometer with 
photon-number-resolving coincidence detection. In a re- 
cent paper we reported on the analysis of a coherent 
state using a single photon-number-resolving detector 
[23( | . We have extended this experimental capability to 
two simultaneously operating visible light photon coun- 
ters (VLPCs) [24j, cryogenic photodetectors that provide 
a current pulse of approximately 40,000 electrons per de- 
tected photon. The VLPCs were maintained at 8 K in a 
helium flow cryostat, and their photocurrent was ampli- 
fied by low- noise, room temperature amplifiers. We mea- 
sured a detection efficiency of 35% and a dark count rate 
of 3 x 10 5 for each detector under our operating condi- 
tions. Custom electronics processed the amplified VLPC 
current pulses to perform gated, fast coincidence detec- 
tion. We were thus able to determine, for each pulse, how 
many photons were detected at both ports c and d. A 
Ti:sapphire pulsed laser, attenuated such that n = 1.08 
photons, provided the input state. Since a coherent state 
maintains its form under linear loss [ll|, the presence of 
loss after the interferometer is completely equivalent to 
a lossless interferometer fed by a weaker input state. We 
use ft to signify the average number of photons in the 
detected state per pulse, after all losses. We were lim- 
ited by the amplifiers to measuring up to four photons 
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FIG. 4: (color online). Difference between the mean value 
of the phase estimator, (O es t), obtained after 150 replica of 
p = 1000 independent measurements and the true value of the 
phase shift 6. The vertical bars are the mean square fluctua- 
tions al st = ((Best - (Best)) 2 )- The result a es t > (6 - Best) 
proves that our protocol provides an unbiased experimental 
phase estimation. 



per pulse Fig. Hfb), but at n = 1.08 the probability of 
detecting five or more photons is negligible. The phase 
shift 9 was changed by tilting a birefringent crystal inside 
the interferometer. 

The first part of the experiment consists of the cali- 
bration of the interferometer. At different, known, values 
of the phase shift, we measured N c and Nd for each of 
200, 000 laser pulses. This procedure allows us to deter- 
mine experimentally both P{N C , Nd\4>) and P((j)\N Cl Nd). 
In Fig. ([2]) we compare the ideal and the experimental 
phase distributions. The agreement is quite good, and 
the discrepancies can be attributed to imperfect photon- 
number discrimination by the detectors. To fit the data 
in Fig.©, we introduce the probability P(N C , N d \N' c , N' d ) 
to measure N c and Nd when N' c and N' d photons were re- 
ally present. Taking this into account, the experimental 
phase probability distribution is 



PMN^Nd) 



N'N' 



\N' c ,N f d )P{N f c ,N f d \N c ,N d ), 



(5) 

where P((j>\N' c , N d ) are the ideal probabilities Eq.gj). The 
weights P(N' c1 N' d \N Cl Nd) can be retrieved from a fit of 
the experimental calibration distributions P(0|iV c , N d ), 
see Fig.®. The quantity P(N C , N d \N c , N d ), equal to one 
in the ideal case, is 0.54 in Fig.([2]A) (corresponding to 
the worst case among all distributions), 0.67 in ([2j3) and 
(00), and 0.87 in {2p). 

After the calibration, we can proceed with the 
Bayesian phase estimation experiment. For a certain 
value of the phase shift, we input one laser pulse and 
detect the number of photons N c and N d . We repeat 
this procedure p times obtaining a sequence of indepen- 
dent results {Nc Z \N^}i=i m „p. The p photon-number 
measurements comprise a single phase estimation. The 
overall phase probability is given by the product of the 



FIG. 5: (color online). Phase sensitivity as a function of 
the true value of the phase shift. Circles are obtained from 
Bayesian distributions Eq.JS}, with p = 1000. The error bars 
give the fluctuations of y/pAO obtained with 150 indepen- 
dent replica of the experiment. The solid black line is the 
theoretical prediction Eq.j3}. The dashed blue line is the 
CRLB calculated with the experimental distributions. Tri- 
angles are the uncertainty obtained with a generalization of 
the estimator Eq.(pQ) taking into account the experimental im- 
perfections, while the dotted red line is the phase sensitivity 
predicted by Eq.©. 



distributions associated with each experimental result: 



P{4>\{N^,N^} i=1 ... p ) oc f[P Rt (<J>\N?\N 
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The phase estimator est is given by the mean value of 
the distribution, 6 est = dt/x/) P ((/>\{N^\ A^ } } i=1 ... p ) , 
and the phase uncertainty A6 is the 68.27% con- 
fidence interval around O es t- An example of 
P(0|{JVc (<) ,JVi O }»=i...p) is given in Fig.©, for 0/tt = 0.24 
and for different values of p. Since the average number of 
photons of the coherent input state is small, for p ~ 1 the 
phase uncertainty is of the order of the prior knowledge, 
AO 7r. As p increases, the probability distribution be- 
comes Gaussian and the sensitivity scales as A0 oc 1/ ^/p, 
in agreement with the central limit theorem. 

In Fig.(j4j), we show the difference between the mean 
value of the phase estimator (0 es t)> obtained from 150 
phase estimations, each with p = 1000, and the true value 
of the phase shift. The bars are the mean square fluctua- 
tion. The important result is that our protocol provides 
an experimentally unbiased phase estimation over the en- 
tire phase interval. 

The main result of this Letter is presented in Fig. ([5]). 
We show the phase sensitivity for different values of the 
phase shift calculated from the distribution Eq.© 
with p = 1000 photon- number measurements. The 
circles are the mean value of A0, and the bars give 
the corresponding mean square fluctuation, obtained 
from 150 independent phase measurements. The dashed 
blue line is the CRLB calculated with the experimen- 
tal probability distributions, A#fi t 
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0.1 < 6/tt < 0.9, it follows the theoretical prediction 
(solid black line), Eq.(j3]), where n = 1.08 has been in- 
dependently calculated from the collected data. Around 
= 0, 7r, where the photons have higher probability to 
exit through the same output port, AOm increases as a 
consequence of the decreased sensitivity of our detectors 
to higher photon number states. Even though the phase 
sensitivity of our apparatus becomes worse near = 0, 7r, 
it never diverges. Triangles show the phase uncertainty 
obtained with Eq.([Tj), but taking into account the exper- 
imental noise. The estimator is obtained by inverting 
the equation M p = ncos(a + 6) + b {2^. This strategy 
provides an unbiased estimation with a sensitivity close 
to the one predicted by Eq.([2|) (dotted red line). The 
reason for the superior performance of the Bayesian pro- 
tocol can be understood by noticing that, in Eq.([2|), the 
phase estimate is retrieved only from the measurement 
of the photon number difference, which, simply does not 
exploit all of the available information. 

In [lj] YMK first proposed a generalization of the 
estimator Eq.(pp) to take into account the whole infor- 
mation in the output measurements. Their estimator, 
©ymk = arccos[(7V c — Nd)/(N C + A^)], gives a phase in- 
dependent sensitivity, A0 = l/^/N c + N^. Notice that 
@ymk coincides with the Maximum Likelihood estima- 
tor in the ideal, noiseless, MZ interferometer. However, 
it is not obvious how to generalize Oymk for real inter- 
ferometry, where classical noise is present and the YMK 
estimator is different from Oml- In general, because of 
correlations between N c — and N c + this estima- 
tor becomes strongly biased in the presence of noise as 
we have verified using our experimental data [27[ . More- 
over, the YMK estimator cannot be extended when both 
input ports of the interferometer are illuminated. Con- 
versely, the Bayesian analysis holds for general inputs 
and, in particular, it predicts a phase-independent sen- 
sitivity when squeezed vacuum is injected in the unused 
port of the MZ [27[ , which reaches a sub shot-noise sensi- 
tivity [13]. It should be noted that detection losses again 
become important when attempting to use nonclassical 
light to overcome the shot-noise limit Eq.(j3]). 

In [jjl, Hradil et al. used a Bayesian approach for 
a Michelson-Morley neutron interferometer (single out- 
put detection) and discussed theoretically the MZ. Their 
analysis was based on specific assumptions about the in- 
terferometric classical noise which are not satisfied in the 
case discussed in this Letter. Different approaches with 
adaptive measurements [20| and positive operator value 
measurements [To| have been also suggested. While these 
strategies might be important for interferometry at the 
Heisenberg limit, they are not necessary in our case. 

In conclusion, we have presented a Bayesian phase es- 
timation protocol for a MZ interferometer fed by a single 
coherent state. The protocol is unbiased and provides 
a phase sensitivity that saturates the ultimate Cramer- 
Rao uncertainty bound imposed by quantum fluctua- 
tions. We have been able to implement the protocol 
with two photon-number-resolving detectors at the out- 



put ports of a weakly-illuminated interferometer. Yet, 
the method can be generalized to the case of high inten- 
sity laser interferometry and photodiode detectors. In 
this case, the limit Eq.(|3]) becomes harder to achieve be- 
cause of larger electronic noise and lower photon number 
resolution, however it should still be possible to demon- 
strate a phase independent sensitivity. Our results are of 
importance to quantum inference theory and show that 
the MZ interferometer does not require phase-locking in 
order to reach an optimal sensitivity. 

Acknowledgment. This work has been partially sup- 
ported by the US DoE and NSF Grant No. 0304678. 



[1] L. Mach, Zeitschr. f. Instrkde. 12, 89 (1892). 

[2] L. Zehnder, Zeitschr. f. Instrkde. 11, 275 (1891). 

[3] Y. Ji, et al, Nature 422, 415 (2003). 

[4] Atom Interferometry. P.R. Berman edt., Academic Press, 
New York, 1997. 

[5] Y. Wang et al, Phys. Rev. Lett. 94, 090405 (2005). 

[6] F. Zenhausern, Y. Martin & H.K. Wickramasinghe, Sci- 
ence 269, 1083 (1995). 

[7] A. Peters, K Y. Chung & S. Chu, Nature 400, 849 (1999). 

[8] M.J. Holland & K. Burnett. Phys. Rev. Lett. 71, 1355 
(1993). 

[9] L. Pezze & A. Smerzi, Phys. Rev. A 73, 011801(R) 
(2006). 

[10] B.C. Sanders & G.J. Milburn, Phys. Rev. Lett. 75, 2944 
(1995). 

[11] M.O. Scully & M.S. Zubairy, Quantum Optics, Cam- 
bridge University Press 1997. 

[12] J.P. Dowling, Phys. Rev. A 57, 4736 (1998). 

[13] CM. Caves, Phys. Rev. D 23, 1693 (1981). 

[14] B. Yurke, S.L. McCall & J.R. Klauder, Phys. Rev. A 33, 
4033 (1986). 

[15] H. Cramer. Mathematical methods of statistics, Prince- 
ton, Princeton university press, 1946. 

[16] C.R. Rao, Bull. Calcutta Math. Soc. 37, 81 (1945). 

[17] R.A. Fisher, Proc. Camb. Phi. Soc. 22, 700 (1925). 

[18] C.W. Helstrom, Quantum Detection and Estimation 
Theory Academic Press, New York, 1976. 

[19] Z. Hradil, et al, Phys. Rev. Lett. 76, 4295-4298 (1996). 

[20] D.W. Berry & H.M. Wiseman, Phys. Rev. Lett. 85, 5098 
(2000). 

[21] J.W. Noh, A. Fougeres & L. Mandel, Phys. Rev. Lett. 
67, 1426 (1991); J. Rehacek, et al, Phys. Rev. A 60, 473 
(1991). 

[22] The Quantum Mechanical probability to detect N c and 
Nd particles at the output port c and d, respectively, 
is given by P(N c ,N d \ct>) = ^(N^N^e^y |Vw}| 2 , 
where J y = (a^b — ¥a)/2i, with a, b bosonic annihilation 
operators. For the configuration studied in this Letter 

W. N dW = J^ggf ( - ± ) 2JV < ( cos f . 

[23] G. Khoury, et al, Phys. Rev. Lett. 96, 203601 (2006). 

[24] G.B. Turner, et al, Proceedings of the Workshop on Scin- 
tillating Fiber Detectors, Notre Dame University, edited 
by R. Ruchti (World Scientific, Singapore, 1994). 

[25] A fit of the calibration data gives (N c —Nd) (0) = ft cos(#+ 
a) + b, with a = 0.024 and b = 0.008. 

[26] The distributions P fit (iVi, N 2 \0) are obtained from Eq.© 
with the Bayes theorem. 

[27] Manuscript in preparation. 



